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ABSTRACT 

Aims.  We  aim  at  measuring  mass-loss  rates  and  the  luminosities  of  a  statistically  large  sample  of  Galactic  bulge  stars  at  several 
galactocentric  radii.  The  sensitivity  of  previous  infrared  surveys  of  the  bulge  has  been  rather  limited,  thus  fundamental  questions 
for  late  stellar  evolution,  such  as  the  stage  at  which  substantial  mass-loss  begins  on  the  red  giant  branch  and  its  dependence  on 
fundamental  stellar  properties,  remain  unanswered.  We  aim  at  providing  evidence  and  answers  to  these  questions. 

Methods.  To  this  end,  we  observed  seven  15x15  arcmin2  fields  in  the  nuclear  bulge  and  its  vicinity  with  unprecedented  sensitivity 
using  the  IRAC  and  MIPS  imaging  instruments  on-board  the  Spitzer  Space  Telescope.  In  each  of  the  fields,  tens  of  thousands  of  point 
sources  were  detected. 

Results.  In  the  first  paper  based  on  this  data  set,  we  present  the  observations,  data  reduction,  the  final  catalogue  of  sources,  and  a 
detailed  comparison  to  previous  mid-IR  surveys  of  the  Galactic  bulge,  as  well  as  to  theoretical  isochrones.  We  find  in  general  good 
agreement  with  other  surveys  and  the  isochrones,  supporting  the  high  quality  of  our  catalogue. 

Key  words.  Galaxy:  bulge  —  Galaxy:  stellar  content  —  Infrared:  stars  —  Stars:  late-type  —  Stars:  mass-loss  —  Stars:  AGB  and 
post-AGB 


1.  Introduction 


The  Galactic  bulge  (GB),  an  important  dynamical  and  morpho¬ 
logical  component  of  our  Galaxy,  offers  an  environment  distinct 
from  the  Galactic  disk  for  study  of  stellar  populations,  stellar 
evolution,  and  the  mass-loss  processes  that  accompany  and,  in 
the  end,  control  the  last.  Understanding  and  calibrating  the  phys¬ 
ical  processes  whereby  mass  ejected  by  evolved  stars  into  the 
bulge  environment  is  recycled  back  into  new  generations  of  stars 
requires  a  statistical  knowledge  of  mass  loss  as  a  function  of  fun¬ 
damental  stellar  parameters  in  this  region.  Because  of  the  limited 
sensitivity  of  previous  surveys  of  the  bulge,  fundamental  ques¬ 
tions  for  late  stellar  evolution,  such  as  the  stage  at  which  substan¬ 
tial  mass-loss  begins  on  the  red  giant  branch  (RGB),  and  its  de¬ 
pendence  on  fundamental  stellar  properties,  remain  unanswered. 
The  GB  is  an  ideal  laboratory  for  addressing  these  issues,  provid¬ 
ing  a  very  large  sample  of  stars  at  an  almost  identical  distance. 

We  therefore  observed  seven  15  x  15  arcmin2  fields  that  sam¬ 
ple  a  range  of  distances  from  the  Galactic  centre  with  unprece¬ 
dented  sensitivity  using  the  Infrared  Array  Camera  (IRAC;  Fazio 
et  al.  |2004|>  and  the  Multiband  Imaging  Photometer  for  Spitzer 


*  Table  5  as  well  as  the  fits  mosaics  are  only  available  in  elec¬ 
tronic  form  at  the  CDS  via  anonymous  ftp  to  cdsarc.u-strasbg.fr 
(130.79.128.5)  or  via  http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 


(MIPS;  Rieke  et  al.|2004|l,  the  imaging  instruments  on-board  the 
Spitzer  Space  Telescope  (Werner  et  al.  |2004|l,  in  order  to  de¬ 
termine  mass-loss  rates  and  luminosities  of  a  statistically  large 
sample  of  stars  at  several  galactocentric  radii.  These  data  en¬ 
able  us  to  detect  stars  with  very  low  mass-loss  rates  through 
their  infrared  excess,  determine  the  dependence  of  the  mass- 
loss  rate  on  luminosity  and  effective  temperature  along  the  giant 
branches,  and  conduct  a  census  of  mass-losing  stars  at  different 
rates.  The  observations,  together  with  existing  studies  that  probe 
higher  mass-loss  rate  stars,  will  enable  us  to  infer  the  total  rate 
of  mass  loss  in  the  bulge,  a  key  input  to  evolutionary  models  of 
the  bulge.  The  data  have  already  led  to  the  discovery  of  mid-IR 
log  P  vs.  magnitude  relations  (Glass  et  al.  2009 1.  In  this  paper 
we  present  the  observations,  the  data  reduction,  and  the  source 
catalogue.  We  also  compare  the  data  to  previous  mid-IR  surveys 
of  the  bulge  and  to  theoretical  isochrones. 


The  outline  of  this  paper  is  as  follows.  In  ^2]  we  describe 
the  observations,  the  data  reduction,  the  different  steps  of  point 
source  extraction,  and  how  the  catalogues  were  created.  We  then 
proceed  in  ^3]  with  checking  our  photometric  data  against  other 
mid-IR  surveys  of  the  Galactic  bulge.  Further  checks  are  pre¬ 
sented  in  ^4]  where  colour-magnitude  diagrams  (CMDs)  are 
compared  to  theoretical  isochrones.  Finally,  ^5]  draws  conclu¬ 
sions  on  the  catalogues  and  the  data  quality. 
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2.  Observations  and  Data  Reduction 

2.1.  Field  selection  and  observations 


The  locations  of  the  observed  fields  were  chosen  to  sample  the 
bulge  on  a  variety  of  scales,  to  measure  how  the  mix  of  stellar 
populations  varies  with  Galactic  latitude.  They  were  also  chosen 
to  avoid  the  relative  intense,  saturating  emission  from  near  the 
Galactic  plane.  Our  innermost  fields  are  within  the  central  stel¬ 
lar  cusp,  which  presumably  contains  stars  of  various  ages  (c.f. 
|Blum  et  al.|2003[|Figer  et  akp004).  Because  the  stars  within  this 
domain  are  believed  to  have  formed  within  the  central  molecu¬ 
lar  zone  and  then  to  have  diffused  into  an  increasingly  thicker 
distribution  as  a  result  of  scattering  off  molecular  clouds  (|Kim| 
&  MorrispOOl  i,  we  expect  a  vertical  segregation  of  stellar  ages. 


Thus,  inasmuch  as  the  luminosities  and  mass-loss  rates  of  red 
giants  and  asymptotic  giant  branch  (AGB)  stars  depend  on  their 
masses,  hence  their  ages,  the  radial  distributions  of  the  different 
kinds  of  evolved  stellar  objects  can  be  used  to  model  the  star 
formation  and  dynamical  history  of  this  region. 

Two  fields  were  selected  to  sample  the  nuclear  bulge  at 
(lb)  =  ((TOO,  -LOO)  and  (l,b)  =  (0163,-0136),  namely  N1 
and  N  2,  both  of  which  were  observed  by  ISOGAL  (project 
for  imaging  part  of  the  Galaxy  using  ISO,  the  Infrared  Space 
Observatory;  Omont  et  al.  2003j>  at  7  and  1 5  gm.  Five  fields 
were  selected  beyond  the  nuclear  bulge  in  areas  where  they  over¬ 
lap  the  Optical  Gravitational  Lensing  Experiment  III  (OGLE-IH; 
tJdalski  et  al.|2000|>  micro-lensing  survey.  Four  of  these  are  lo¬ 


cated  below  the  Galactic  plane  along  a  radial  vector  that  sub¬ 
tends  the  minor  axis  of  the  Galaxy  at  an  angle  of  about  1315, 
and  terminates  at  the  well-studied  field  in  Baade’s  window  at 
(l,  b)  -  ( 1703,  -3183),  containing  the  globular  cluster  NGC  6522 
(e.g.  Glass  et  al."fl999)l.  The  three  fields  inside  this  are  located 
at  (l,b)  =  (0130,  -1142),  (0156,  -2723),  and  (0176,  -3107).  The 
fifth  field  is  positioned  above  the  plane  at  (lb)  =  (2187,0135). 
All  fields  are  approximately  rectangular  in  right  ascension  (RA) 
and  declination  (Dec)  coordinates.  Figure  [T| shows  the  location 
of  the  seven  observed  fields  with  respect  to  the  Galactic  cen¬ 
tre,  and  Table  |T] summarises  some  of  their  main  characteristics. 
The  range  in  RA  and  Dec  (J2000)  given  in  Table  |T]  refers  to 
where  there  is  full  overlap  between  MIPS  and  all  four  IRAC 
bands.  Particularly  in  IRAC  we  do  have  some  coverage  outside 
the  given  range. 

Observations  of  the  GB  fields  were  performed  using  the 
IRAC  instrument  in  all  four  channels  at  3.6, 4.5,  5.8,  and  8.0 pm, 
and  the  MIPS  instrument  in  the  24  pm  channel  on  board  the 
Spitzer  Space  Telescope  within  the  programme  2345.  Bulge  field 
N  2  was  observed  only  with  MIPS  in  our  programme,  because 
this  field  was  covered  by  IRAC  observations  within  the  gen¬ 
eral  observer  (GO)  programme  3677  (principal  investigator:  S. 
Stolovy;  Ramirez  et  al.|2008|.  The  IRAC  observations  were  car¬ 
ried  out  on  2005  March  30  and  31,  and  the  MIPS  observations 
between  2005  April  08  and  13. 


With  IRAC,  observations  were  done  in  the  full-array  read 
out-mode,  one  frame  per  pointing  with  2  s  integration  time  per 
frame  -  the  shortest  full-array  integration  time  was  used  to  min¬ 
imise  the  effects  of  saturation  and  resulting  latency  problems. 
The  mapping  was  done  in  a  6  x  6  rectangular  grid  with  step  size 
260",  with  five  dither  positions,  and  a  medium  scale  factor,  giv¬ 
ing  a  total  exposure  time  of  10  s  per  pixel.  The  MIPS  observa¬ 
tions  were  obtained  using  the  photometry  /raster  mode  with  10  s 
integration  time  and  full-array  read-out  mode,  3x3  rectangular 
grid  and  two  cycles,  giving  a  total  integration  time  of  331  s  per 
pixel. 
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Fig.  1.  Map  of  the  seven  fields  observed  in  this  campaign  with 
the  labels  used  throughout  this  paper. 


The  fields  of  view  of  IRAC  channels  1  and  3,  2  and  4,  as  well 
as  MIPS  did  not  fully  overlap.  Nevertheless,  the  MIPS  field  of 
view  is  fully  contained  in  all  four  IRAC  fields  of  view. 


2.2.  Data  reduction 


We  corrected  the  IRAC  basic  calibrated  data  (bed)  files  to  mit¬ 
igate  artifacts  such  as  muxbleed  or  column  pulldown  with  the 
tools  provided  by  S.  Carey  at  the  Spitzer  Science  Center  (SSC). 
Post-bed  processing  was  then  conducted  on  the  corrected  bed 
files  using  the  MOsaicker  and  Point  source  Extractor  (MOPEX) 


software  and  its  subsystem  APEX  (version  18.2.0,  Markovoz  & 
|Marleau|2005]l. 


2.2.1 .  Step  1 :  Point  source  detection  and  extraction 

Before  further  processing,  mosaics  were  created  by  MOPEX 
for  each  field  and  channel  from  the  corrected  bed  frames.  In 
this  step,  the  pipeline  interpolates  the  input  images  onto  the 
output  grid,  taking  geometric  distortion  into  account.  An  out¬ 
lier  detection  scheme  flags  bad  pixels  and  any  pixels  affected 
by  cosmic  ray  hits  or  moving  objects,  and  these  pixels  are  re¬ 
computed  before  co-addition.  Finally,  the  interpolated  images 
are  co-added  to  one  mosaic  image.  After  mosaicking,  APEX  de¬ 
termines  the  background  by  calculating  the  median  in  a  45  x  45 
pixel  box  around  each  pixel,  and  subtracts  it  from  the  image. 
These  background-subtracted  images  were  used  in  the  detec¬ 
tion  step.  Then,  background  fluctuations  in  the  images  were  esti¬ 
mated  and  noise  images  derived,  which  were  used  for  signal-to- 
noise  ratio  (SNR)  estimation  and  the  generation  of  point  source 
probability  images.  With  those,  the  detection  table  was  com¬ 
piled.  We  chose  a  detection  threshold  of  2.3cr  above  the  back¬ 
ground.  We  used  both  the  point  response  function  (PRF)  fitting 
capability  and  the  aperture  photometry  functionality  with  circu¬ 
lar  apertures  with  radii  of  2,  3,  and  4  pixels.  Larger  apertures 
were  not  applicable  owing  to  the  crowding  in  the  fields.  The  tiles 
used  for  PRF  fitting  also  had  to  be  chosen  to  be  very  small  (3x3 
pixels).  We  used  the  most  up-to-date  mean  PRFs  as  provided  by 
the  SSC,  and  set  the  PRF  normalisation  radius  accordingly. 

The  measured  x~  values  in  our  data  seem  to  be  very  high, 
even  for  successfully  fitted  sources,  e.g.  in  the  least-crowded 
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Table  1.  Main  characteristics  of  the  seven  fields  towards  the  GB  observed  with  Spitzer  as  part  of  programme  2345. 


Field  name 

RA  range  (h  m  s) 

Dec  range  (°  '  ") 

1  centre 

b  centre 

Date  IRAC  obs 

Date  MIPS  obs 

Bulge  1 

17  51  14.2 

..17  52  36.7 

-29  15  22 

..-29  34  15 

0.30 

-1.42 

2005  03  30 

2005  04  09 

Bulge  2 

17  55  02.1 

. .  17  56  24.2 

-29  26  25 

. .  -29  45  13 

0.56 

-2.23 

2005  03  31 

2005  04  10 

Bulge  3 

17  58  50.7 

..18  00  13.3 

-29  41  08 

..-29  59  59 

0.76 

-3.07 

2005  03  30 

2005  04  10 

Bulge  4 

17  50  16.1 

..17  51  36.4 

-26  08  56 

..-26  27  51 

2.87 

0.35 

2005  03  30 

2005  04  09 

Bulge  N  1 

17  48  51.9 

..17  50  14.8 

-29  17  42 

. .  -29  36  42 

0.00 

-1.00 

2005  03  30 

2005  04  08 

Bulge  N  2 

17  47  55.8 

..17  49  16.3 

-28  28  41 

..-28  46  19 

0.63 

-0.36 

- 

2005  04  13 

NGC  6522 

18  02  30.4 

..18  03  53.5 

-29  49  26 

. .  -30  08  24 

1.03 

-3.83 

2005  03  30 

2005  04  10 

field  NGC  6522,  they  are  between  2.5  and  3.3  in  the  IRAC  chan¬ 
nels  1  and  2,  between  1 .4  and  1 .8  in  IRAC  3,  and  only  in  IRAC  4 
at  the  perfect  value  of  1 .  This  may  stem  from  confusion  noise, 
which  is  not  included  in  the  provided  uncertainty  images  and/or 
from  the  high  source  density.  Confusion  noise  is  created  by 
the  amplitude  variations  from  PRFs  of  closely  spaced  sources 
( Hacking  &  Houck|  1 987[  Rieke  et  al.|1995 1.  Modelling  the  con¬ 
fusion  noise  has  not  improved  the  resulting  x2  values.  This  has 
consequences  for  the  reliability  of  automatic  de-blending  proce¬ 
dures,  which  were  not  working  properly  in  this  field.  We  thus 
had  to  disable  the  active  and  passive  de-blending  capabilities  of 
APEX  (see  Appendix  [A}. 


In  all  MIPS  fields,  we  first  detected  bright  sources,  removed 
the  Airy  rings  of  these  sources  to  do  the  source  detection  on  the 
residual  images,  and  performed  the  point  source  extraction  on 
the  original  images.  The  southern  half  of  the  field  N  2  has  strong 
diffuse  emission.  In  that  area,  the  detection  threshold  is  certainly 
higher  than  in  areas  with  low  diffuse  background  emission,  as  is 
also  evident  from  the  histograms  of  magnitudes  in  Appendix [B| 


2.2.2.  Step  2:  Assigning  quality  grades  to  the  sources 

Three  different  effects  can  be  seen  when  plotting  the  flux  ob¬ 
tained  by  PRF  fitting  versus  that  found  with  aperture  photome¬ 
try  for  all  the  extracted  sources  (Fig.|2|.  The  first  effect  is  present 
when  we  use  a  large  aperture  with  a  radius  of  3  or  4  pixels.  At 
very  low  PRF  fluxes,  there  are  several  sources  that  show  dispro¬ 
portionately  large  aperture  fluxes,  i.e.  are  above  the  dashed  line. 
These  sources  are  faint  sources  close  to  a  bright  source.  In  case 
of  the  larger  aperture,  we  therefore  already  get  a  significant  con¬ 
tribution  in  the  aperture  from  the  bright  source.  The  second  effect 
is  relevant  for  some  bright  sources,  whose  PRF  fluxes  are  smaller 
than  the  aperture  fluxes.  These  sources  are  super-saturated;  i.e., 
the  radial  intensity  profile  has  a  dip  in  the  centre  and  so  it  shows 
a  double  peak  structure.  In  a  few  cases,  depending  on  the  depth 
of  the  central  valley,  these  two  peaks  are  identified  as  two  indi¬ 
vidual  sources  that  are  each  fitted  with  a  fainter  PRF;  i.e.,  both 
sources  lie  above  the  dashed  line.  In  other  cases,  the  fitted  PRF 
ends  up  with  a  central  flux  between  that  of  the  central  valley  and 
that  of  the  wings,  meaning  PRF  fitting  underestimates  the  total 
flux,  so  these  sources  are  also  shifted  to  the  left  in  Fig.  [2] 

The  SSC  homepage  provides  values  of  the  maximum  flux 
of  unsaturated  point  sources  as  a  function  of  integration  time. 
Table[2]presents  the  flux  values  interpolated  to  a  frame  exposure 
time  of  1 .2  s,  along  with  a  lower  flux  limit.  This  lower  flux  limit 
was  adopted  on  subjective  grounds  for  the  quality  grade  labelling 
(see  below).  We  thus  have  to  be  cautious  with  sources  brighter 
than  about  320  mJy  in  IRAC  1 .  Consistent  with  this  warning  is 
a  third  slight  effect  that  seems  to  set  in  at  fluxes  close  to  this 
value  and  affects  all  brighter  sources:  the  PRF  fluxes  increase 
disproportionately  compared  to  the  aperture  fluxes,  leading  to  a 
bending  to  the  right,  away  from  the  dashed  line  in  Fig.  [2]  (top). 


Table  2.  Flux  limits  applied  for  assigning  quality  grades. 


Band 

Lower  limit  (mJy) 

Saturation  limit  (mJy) 

IRAC  3.6yum 

40 

320 

IRAC  4.5  pm 

30 

330 

IRAC  5.8/rm 

20 

2300 

IRAC  8.0/rm 

20 

1200 

MIPS  24  pm 

40 

220 

Aperture  photometry  thus  measures  the  saturated  plateau  of  the 
PRF,  while  PRF  fitting  seems  to  ignore  the  missing  flux  that 
would  be  present  beyond  the  plateau  and  fits  the  wings.  The  indi¬ 
vidual  source  then  lies  below  the  dashed  line,  with  a  larger  PRF 
flux  than  the  aperture  flux. 

If  dividing  the  aperture  flux  by  the  linear  least  squares  fit  to 
sources  with  x 2  <  5,  the  fluxes  measured  with  both  ways  are 
almost  identical  for  many  of  the  sources  with  x2  <  5.  The  slope 
of  this  fit  is  not  the  aperture  correction  factor,  to  be  derived  in 
Sect.  2.2.3  Except  for  the  sources  affected  by  the  three  afore¬ 
mentioned  effects,  the  relative  difference  calculated  as 


f  ^corrected 

IPRF  ^aDerture 

rel.  difference  = - - -  (1) 

fpRF 

scatters  around  7%  (Fig. [2]  bottom).  The  superscript  “corrected” 
in  Eq.[T]  indicates  that  the  aperture  flux  has  been  divided  by  the 
linear  fit  to  sources  with  x2  <  5.  Sources  that  are  within  this 
range  are  considered  as  quality  grade  A.  This  test  is  performed 
first,  therefore  grade  A  sources  span  the  whole  range  of  fluxes, 
from  faint  to  almost  saturated  sources. 

Faint  sources  with  fluxes  smaller  than  a  certain  lower  limit 
(Table  [2}  and  outside  the  7%  range  are  classified  as  quality 
grade  B.  Sources  whose  fluxes  are  above  the  saturation  limit  are 
classified  as  quality  grade  C.  They  would  have  to  be  treated  dif¬ 
ferently  with  dedicated  tools.  However,  since  these  sources  are 
probably  mostly  foreground  stars,  they  are  not  of  prime  interest 
to  the  science  goals,  thus  we  decided  not  to  invest  more  time 
in  their  flux  determination.  The  interested  reader  is  invited  to 
download  and  analyse  the  mosaics,  which  will  be  made  avail¬ 
able  as  on-line  material  at  Centre  de  Donnees  astronomiques  de 
Strasbourg  (CDS). 

Sources  that  failed  all  of  the  above  criteria  are  classified  as 
quality  grade  D.  These  are  blended  sources  for  which  the  PRFs 
overlap  such  that  even  the  aperture  with  the  smallest  radius  in¬ 
cludes  both  peaks;  i.e.,  no  reliable  flux  measurement  is  available 
for  these  sources.  This  effect  is  the  same  as  the  first  one  men¬ 
tioned  above  related  to  the  larger  aperture  with  a  radius  of  3  pix¬ 
els  or  larger. 

The  statistics  of  grades  for  all  seven  fields  and  all  five  bands 
are  given  in  Table  [3]  The  percentage  given  for  grades  A  and  B 
is  relative  to  the  total  number  of  sources  in  the  catalogue  for  the 
respective  field  that  is  given  in  the  last  column.  The  seventh  col¬ 
umn  gives  the  number  of  sources  not  detected  in  the  respective 
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Table  3.  Statistics  of  sources  with  their  grades  for  all  seven  fields  and  five  bands. 


Field 

Band 

Grade  A 

Grade  B 

Grade  C 

Grade  D 

Not  detected 

total 

Bulge  1 

3.6pm 

12446(18%) 

39177  (57%) 

183 

626 

16661 

69093 

4.5  gm 

16173  (23%) 

35598  (52%) 

85 

216 

17021 

5.8pm 

13839  (20%) 

21757  (31%) 

1 

169 

33327 

8.0pm 

12669(18%) 

17044  (25%) 

7 

96 

39277 

24  pm 

577  (  1%) 

1197(2%) 

1 

2 

67316 

Bulge  2 

3.6pm 

14889  (25%) 

30525  (51%) 

197 

748 

13963 

60322 

4.5  pm 

17495  (29%) 

28190  (47%) 

72 

492 

14073 

5.8pm 

13712(23%) 

17900  (30%) 

0 

214 

28496 

8.0pm 

11960  (20%) 

15752  (26%) 

1 

152 

32457 

24  pm 

503  (  1%) 

1119(2%) 

4 

4 

58692 

Bulge  3 

3.6pm 

13476  (26%) 

27927  (53%) 

131 

101 

11031 

52666 

4.5  pm 

16329  (31%) 

25159(48%) 

62 

64 

11052 

5.8pm 

12085  (23%) 

15880  (30%) 

0 

35 

24666 

8.0pm 

10712  (20%) 

13220  (25%) 

4 

24 

28706 

24  pm 

383  (  1%) 

976 ( 2%) 

1 

2 

51304 

Bulge  4 

3.6pm 

12938  (22%) 

33046  (57%) 

136 

197 

12091 

58408 

4.5  pm 

15689  (27%) 

29961  (51%) 

75 

95 

12588 

5.8pm 

11311  (19%) 

18283  (31%) 

3 

177 

28634 

8.0  pm 

9239  (16%) 

13588  (23%) 

19 

98 

35464 

24  pm 

267  ( 0%) 

817  (  1%) 

3 

4 

57317 

Bulge  N  1 

3.6pm 

12553  (19%) 

35826  (55%) 

198 

724 

15534 

64835 

4.5  pm 

15402  (24%) 

34527  (53%) 

85 

296 

14525 

5.8pm 

11819(18%) 

21845  (34%) 

5 

520 

30646 

8.0pm 

12035  (19%) 

16565  (26%) 

11 

144 

36080 

24  pm 

467  (  1%) 

1126(2%) 

3 

23 

63216 

Bulge  N  2 

24  pm 

21  ( 5%) 

316(74%) 

17 

74 

0 

428 

NGC  6522 

3.6pm 

12430  (26%) 

24924  (53%) 

79 

55 

9796 

47284 

4.5  pm 

15637  (33%) 

21263  (45%) 

32 

21 

10331 

5.8pm 

12011  (25%) 

11907  (25%) 

0 

30 

23336 

8.0  pm 

10263  (22%) 

10432  (22%) 

0 

12 

26577 

24  pm 

431 (  1%) 

1152(2%) 

0 

0 

45701 

channel,  but  detected  in  at  least  one  of  the  other  bands  plus  in 
the  reference  catalogues  (see  Sect.  2.2.4). 


2.2.3.  Step  3:  Determination  of  the  aperture  corrections 

Since  we  do  not  use  a  larger  aperture  to  estimate  the  background, 
which  would  give  unreliable  photometry  in  the  crowded  fields, 
but  do  aperture  photometry  on  the  background-subtracted  im¬ 
ages  instead,  the  aperture  corrections  are  expected  to  be  different 
from  the  standard  values  given  in  the  instrument  handbooks. 

To  derive  aperture  corrections,  we  created  point  spread 
functions  with  stinytim2.0  (|Kristj[2006]>  for  all  five  bands.  We 
then  inserted  artificial  sources  with  known  flux  in  our  mo¬ 
saicked,  background-subtracted  images  and  compared  the  result¬ 
ing  fluxes  given  by  APEX  with  the  initial  fluxes  with  which  we 
had  created  the  sources.  The  slope  of  a  linear  least  squares  fit  to 
the  initial  vs.  measured  flux  gives  the  aperture  correction  factors. 
The  derived  aperture  corrections  for  the  aperture  with  two  pixels 
radius  are  given  in  Table  [4]  along  with  their  statistical  error  as 
derived  from  the  linear  fit.  Also  given  in  the  table  are  the  values 
as  recommended  by  the  SSC,  for  a  2  pixel  radius  on  source  and 
a  10  -  20  pixel  radius  background  annulus  in  the  case  of  IRAC, 
and  a  3'.' 5  on  source  aperture  without  any  background  annulus 
in  the  case  of  MIPS.  The  SSC  recommended  values  are  not  di¬ 
rectly  comparable  to  our  aperture  correction  factors,  so  they  are 
only  given  for  illustrative  purposes.  Despite  the  high  value  of 
the  aperture  correction  factor  for  MIPS  24,  we  adopted  the  small 
aperture  with  2  pixels  radius  because  i)  our  observations  are  not 
photon  noise  limited,  thus  a  larger  aperture  would  not  reduce 
the  error  estimate;  ii)  there  is  a  “population”  of  blended  sources 


Table  4.  Aperture  correction  factors 


Band 

Aperture  correction  SSC  recommended 

IRAC  3.6pm 

1.2129±0.0008 

1.205 

IRAC  4.5  pm 

1.2730±0.0005 

1.221 

IRAC  5.8pm 

1.4693±0.0007 

1.363 

IRAC  8.0  pm 

1.6556±0.0006 

1.571 

MIPS  24  pm 

1.5742±0.0007 

2.560 

whose  flux  increases  disproportionally  with  growing  aperture 
because  more  of  the  neighbouring  source’s  flux  is  measured  with 
larger  aperture;  iii)  the  SSC  lists  aperture  corrections  for  an  even 
smaller  aperture,  so  our  value  is  not  an  extreme.  We  also  cor¬ 
rected  the  PRF  fluxes  with  the  photometric  correction  factors  as 
recommended  by  the  SSC. 

Besides  the  aperture  correction  and  the  photometric  correc¬ 
tion  for  the  PRF  fluxes,  we  did  not  apply  any  other  correction 
factors  such  as  the  colour  correction.  According  to  the  IRAC 
data  handbook,  the  colour  correction  factor  for  blackbody  spec¬ 
tra  with  temperatures  between  2  000  and  5  000  K  is  around  1%, 
except  for  IRAC  4,  where  it  can  reach  as  much  as  2.7%  for  the 
hotter  stars.  Since  the  colour  correction  is  small  and  we  do  not 
know  a  priori  the  nature  and  spectrum  of  our  objects,  we  did  not 
apply  these  corrections. 

2.2.4.  Step  4:  Creating  the  point  source  catalogue 

We  finally  decided  to  adopt  only  the  PRF  fluxes  for  the  final 
catalogue  on  the  basis  of  comparison  with  other  catalogues  and 
with  isochrones  in  CMDs.  With  this  choice,  we  find,  on  the  one 


Stefan  Uttenthaler  et  al. :  Galactic  bulge  giants:  probing  stellar  and  galactic  evolution 


5 


?oo[rr' 

600 1- 

■?  600  L- 
E  E 


”  t 
1  100 t- 


7°°|e . . . ."."M . 'i 


Fig.  2.  Comparison  of  the  fluxes  resulting  from  PRF  fitting 
and  aperture  photometry  for  the  IRAC  1  observations  of  the 
NGC  6522  field.  Top:  PRF  fitting  fluxes  and  aperture  photom¬ 
etry  for  the  aperture  with  a  radius  of  2  pixels;  middle:  the  same 
as  top  panel,  but  for  the  aperture  with  a  radius  of  3  pixels;  bot¬ 
tom:  relative  difference  between  the  measured  fluxes  (3  pixels 
aperture  radius)  versus  the  PRF  flux  for  all  sources.  The  dotted 
lines  in  the  bottom  panel  mark  relative  differences  of  ±7%.  In  the 
upper  and  middle  panels,  the  dashed  line  is  a  linear  least  squares 
fit  to  the  sources  with  <  5,  and  the  solid  horizontal  and  verti¬ 
cal  lines  mark  the  lower  and  saturation  limits  of  Table|3used  for 
quality  grade  labelling. 


hand,  slightly  better  agreement  with  other  catalogues,  and  on 
the  other,  somewhat  reduced  scatter  at  the  faint  end  of  CMDs. 
See  Sect.  [3]  and  [4]  for  a  comparison  with  other  catalogues  and 
isochrones. 

The  last  step  is  the  cross-identification,  for  each  field,  of  the 
extraction  tables  of  all  five  bands  among  each  other  and  with 
other  reference  catalogues.  As  reference  catalogues  we  adopt  the 
Deep  Near  Infrared  Survey  (DENIS;  Epchtein  et  al.|19971>,  the 
2  Micron  All  Sky  Survey  (2MASS;  Skrutskie  et  al.|2006),  the 
Midcourse  Space  Experiment  (MSX)  catalogue,  and  the  Infra- 
Red  Astronomical  Satellite  (IRAS)  point  source  catalogue.  We 
applied  the  following  criteria  to  include  a  source  in  the  final  cat¬ 
alogue.  We  identified  individual  sources  by  their  RA  and  Dec 
position  and  allowed  for  an  error  margin  of  1'.'6,  i.e.  slightly 
more  than  one  pixel  size,  for  the  cross-identification  among  our 
five  bands.  Sources  that  are  detected  in  at  least  two  of  our  five 
bands  are  included  in  any  case  in  the  final  catalogue.  A  cross- 
identification  is  made  for  these  sources  with  a  search  radius  of 
370  with  DENIS  and  2MASS,  and  3070  for  identification  with 
an  MSX  or  IRAS  counterpart.  The  information  about  these  coun¬ 
terparts  is  included  in  the  catalogue.  Sources  that  are  detected  in 
only  one  of  our  five  bands  are  only  included  in  the  final  cata¬ 
logue  if  they  also  have  a  counterpart  in  at  least  one  of  the  ref¬ 
erence  catalogues,  with  the  same  search  radii  applied  as  for  the 
sources  with  two  or  more  Spitzer  detections.  This  procedure  was 
followed  to  include  as  many  of  the  sources  in  the  final  catalogue 
that  were  only  detected  in  MIPS  24.  These  might  very  well  be 
real  objects,  albeit  extremely  red  (i.e.  not  detected  with  IRAC). 
Additionally,  by  following  this  procedure,  more  of  the  sources 
that  are  located  in  the  areas  that  are  not  sampled  by  all  five  bands 
will  be  included  in  the  catalogue.  These  areas  are  also  included 
in  the  final  catalogue.  The  least  reliable  sources  in  our  catalogue 
are  those  with  detection  in  only  one  IRAC  band  and  no  MIPS 
detection,  and  with  (i)  either  no  2MASS  and  DENIS  counterpart 
(depending  on  the  field,  between  3.1%  and  6.7%  of  the  sources), 
or  (ii)  a  DENIS  and  2MASS  counterpart  at  a  distance  between 
1.6  and  3  arcseconds  (depending  on  the  field,  between  0.4%  and 
1.0%  of  the  sources). 

Table  [5]  gives  ten  lines  of  the  band-merged  catalogue  of  the 
Bulge  N  1  field  as  an  example.  The  entire  table  will  be  made 
available  for  download  from  CDS.  A  portion  is  shown  here  for 
guidance  regarding  its  form  and  content. 

The  columns  of  our  point  source  catalogue  are  explained  as 
follows. 


Column  1:  Source  identification,  IAU-conform  identifier. 

Column  2:  Right  ascension  in  degrees  (J2000). 

Column  3:  Declination  in  degrees  (J2000). 

Column  4:  PRF  flux  in  IRAC  1  (in  pjy),  set  to  “-9.99999e+99” 
if  not  detected. 

Column  5:  PRF  flux  in  IRAC  2  (in  pjy),  set  to  “-9.99999e+99” 
if  not  detected. 

Column  6:  PRF  flux  in  IRAC  3  (in  pjy),  set  to  “-9.99999e+99” 
if  not  detected. 

Column  7:  PRF  flux  in  IRAC  4  (in  pjy),  set  to  “-9.99999e+99” 
if  not  detected. 

Column  8:  PRF  flux  in  MIPS  24  (in  /rJy),  set  to 
“-9.99999e+99”  if  not  detected. 

Column  9:  Quality  flags  of  the  IRAC  1,  2,  3,  4,  and  MIPS  24 
channels  (see  Section|2.2.2|),  set  to  if  not  detected. 

Column  10:  Uncertainty  of  the  PRF  flux  in  IRAC  1  (in  piy),  set 
to  “-9.99999e+99”  if  not  detected. 

Column  1 1 :  Uncertainty  of  the  PRF  flux  in  IRAC  2  (in  piy),  set 
to  “-9.99999e+99”  if  not  detected. 
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Column  12:  Uncertainty  of  the  PRF  flux  in  IRAC  3  (in  /rJy),  set 
to  “-9.99999e+99”  if  not  detected. 

Column  13:  Uncertainty  of  the  PRF  flux  in  IRAC  4  (in  g/Jy),  set 
to  “-9.99999e+99”  if  not  detected. 

Column  14:  Uncertainty  of  the  PRF  flux  in  MIPS  24  (in  /.(Jy), 
set  to  “-9.99999e+99”  if  not  detected. 

Column  15:  Observation  flag.  The  first  four  bits  state  if  the  posi¬ 
tion  of  the  source  was  in  the  field  of  view  of  the  IRAC  bands 
(“1”)  or  not  (“0”).  the  fifth  bit  is  for  MIPS  24. 

Column  16:  Distance  to  closest  DENIS  source  (arcsec;  only 

<  370),  set  to  ”  if  no  counterpart. 

Column  17:  DENIS  /-band  magnitude,  set  to  if  no  counter¬ 
part  or  the  counterpart  has  no  /-band  measurement. 

Column  18:  DENIS  /-band  magnitude,  set  to  ”  if  no  counter¬ 
part  or  the  counterpart  has  no  /-band  measurement. 

Column  19:  DENIS  //-band  magnitude,  set  to  if  no  coun¬ 
terpart  or  the  counterpart  has  no  //-band  measurement. 

Column  20:  Distance  to  closest  2MASS  source  (arcsec;  only 

<  370),  set  to  ”  if  no  counterpart.  Only  sources  with  qual¬ 
ity  flags  A,  B,  C,  or  D  in  at  least  one  2MASS  filter  have  been 
considered. 

Column  21:  2MASS  /-band  magnitude,  set  to  if  no  coun¬ 
terpart  within  370  was  found,  or  “-9.999”  if  the  counter¬ 
part  has  no  measured  /-band  magnitude  available  or  too  low 
quality  flag. 

Column  22:  2MASS  //-band  magnitude,  set  to  ”  if  no  coun¬ 
terpart  within  370  was  found,  or  “-9.999”  if  the  counter¬ 
part  has  no  measured  //-band  magnitude  available  or  too  low 
quality  flag. 

Column  23:  2MASS  //-band  magnitude,  set  to  ”  if  no  coun¬ 
terpart  within  370  was  found,  or  “-9.999”  if  the  counter¬ 
part  has  no  measured  //-band  magnitude  available  or  too  low 
quality  flag. 

Column  24:  2MASS  quality  flags,  set  to  if  no  counterpart 
or  too  low  quality  flag. 

Column  25:  Distance  to  closest  IRAS  source  (arcsec;  only 

<  30"),  set  to  ”  if  no  counterpart.  Only  sources  with  qual¬ 
ity  flags  2  or  3  in  at  least  one  IRAS  band  have  been  consid¬ 
ered. 

Column  26:  IRAS  12  /j m  flux  (Jy),  set  to  if  no  counterpart, 
and  to  “-9.9e+99”  if  the  counterpart  has  no  measured  flux 
in  the  12  gun  band  or  too  low  quality  flag. 

Column  27:  IRAS  1 2  gun  quality  flag,  set  to  ”  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  28:  IRAS  25  gm  flux  (Jy),  set  to  if  no  counterpart, 
and  to  “-9.9e+99”  if  the  counterpart  has  no  measured  flux 
in  the  25  gm  band  or  too  low  quality  flag. 

Column  29:  IRAS  25  gm  quality  flag,  set  to  ”  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  30:  IRAS  60  gm  flux  (Jy),  set  to  ”  if  no  counterpart, 
and  to  “-9.9e+99”  if  the  counterpart  has  no  measured  flux 
in  the  60  gm  band  or  too  low  quality  flag. 

Column  3 1 :  IRAS  60  gm  quality  flag,  set  to  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  32:  IRAS  100  gun  flux  (Jy),  set  to  ”  if  no  counterpart, 
and  to  “-9.9e+99”  if  the  counterpart  has  no  measured  flux 
in  the  100 gm  band  or  too  low  quality  flag. 

Column  33:  IRAS  100 gun  quality  flag,  set  to  if  no  counter¬ 
part  or  too  low  quality  flag. 


Column  34:  Distance  to  closest  MSX  source  (arcsec;  only 
<  30"),  set  to  ”  if  no  counterpart.  Only  sources  with  qual¬ 
ity  flags  2,  3,  or  4  in  at  least  one  MSX  band  have  been  con¬ 
sidered. 

Column  35:  MSX  B1  band  flux  (Jy),  set  to  ”  if  no  counter¬ 
part,  and  to  “-9.999e+99”  if  the  counterpart  has  no  mea¬ 
sured  flux  in  the  B1  band  or  too  low  B1  band  quality  flag. 

Column  36:  MSX  B1  band  quality  flag,  set  to  if  no  coun¬ 
terpart  or  too  low  quality  flag. 

Column  37:  MSX  B2  band  flux  (Jy),  set  to  ”  if  no  counter¬ 
part,  and  to  “-9.999e+99”  if  the  counterpart  has  no  mea¬ 
sured  flux  in  the  B2  band  or  too  low  B2  band  quality  flag. 

Column  38:  MSX  B2  band  quality  flag,  set  to  “-”  if  no  coun¬ 
terpart  or  too  low  quality  flag. 

Column  39:  MSX  A  band  flux  (Jy),  set  to  ”  if  no  counterpart, 
and  to  “-9.999e+99”  if  the  counterpart  has  no  measured  flux 
in  the  A  band  or  too  low  A  band  quality  flag. 

Column  40:  MSX  A  band  quality  flag,  set  to  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  41:  MSX  C  band  flux  (Jy),  set  to  ”  if  no  counterpart, 
and  to  “-9.999e+99”  if  the  counterpart  has  no  measured  flux 
in  the  C  band  or  too  low  C  band  quality  flag. 

Column  42:  MSX  C  band  quality  flag,  set  to  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  43:  MSX  D  band  flux  (Jy),  set  to  ”  if  no  counterpart, 
and  to  “-9.999e+99”  if  the  counterpart  has  no  measured  flux 
in  the  D  band  or  too  low  D  band  quality  flag. 

Column  44:  MSX  D  band  quality  flag,  set  to  if  no  counter¬ 
part  or  too  low  quality  flag. 

Column  45:  MSX  E  band  flux  (Jy),  set  to  if  no  counterpart, 
and  to  “-9.999e+99”  if  the  counterpart  has  no  measured  flux 
in  the  E  band  or  too  low  E  band  quality  flag. 

Column  46:  MSX  E  band  quality  flag,  set  to  ”  if  no  counter¬ 
part  or  too  low  quality  flag. 
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3.  Comparison  with  other  missions  and  catalogues 

Our  observations  overlap  with  a  number  of  other  mid-IR  sur¬ 
veys.  These  are 


1. 

2. 


3. 

4. 


ISOGAL:  Survey  with  the  ISO  Camera  (ISOCAM)  on-board 
the  ISO  satellite,  combined  with  DENIS  IJK$  photometry. 
Galactic  Legacy  Infrared  Mid-Plane  Survey  Extraordinaire 
II  (GLIMPSE-II),  a  Spitzer/IRAC  survey  of  the  area  ±1° 
around  the  Galactic  plane  and  the  Galactic  centre  (Benjamin 
let  al.|2003). 


GALactic  CENtre  (GALCEN),  a  Spitzer/IRAC  survey  of  the 
inner  ~  \°A  x  2!0  of  the  Galaxy  (Ramirez  et  al.  2008). 

A  Spitzer/MlPS  survey  of  the  inner  ~  1°5  X  87)  of  the  Galaxy 
(|Hinz  et  al.|2008). 


3.1.  Comparison  with  iSOGAL 


ISOGAL  (Omont  et  al.  2003|>  is  a  survey  with  the  ISOCAM 
instrument  on-board  the  ISO  satellite  in  two  bands  at  7  and 
15  pm,  combined  with  IJKs  photometry  from  the  DENIS  project 
HEpchtein  et  aLl[l997).  The  7  pm  LW2  band  of  ISOCAM  (~ 
5. 0-8. 5  pm,  central  wavelength  6.75  pm)  overlaps  in  wavelength 
with  the  IRAC4  channel  (~  6.4  -  9.4pm,  nominal  wavelength 
7.844pm;  Hora  et  al.  2008).  For  IRAC4,  we  have  overlap  with 
ISOGAL  in  the  fields  Bulge  2,  Bulge  4,  N  1,  and  NGC  6522.  The 
number  of  sources  in  common  (with  less  than  2'.'0  positional 
offset)  is  163  for  Bulge  2,  335  for  Bulge  4,  288  for  N  1,  and 
264  in  NGC  6522.  We  present  here  a  brief  comparison  for  the 
NGC  6522  field,  because  the  ISO  7  pm  photometry  goes  deepest 
in  that  field  (least  crowding  and  diffuse  background  emission). 
The  results  for  the  other  fields  are  very  similar. 

Figure[3]shows  the  magnitude  difference  between  the  Spitzer 
IRAC4  band  and  the  ISO  LW2  band  for  the  NGC  6522  field. 


The  data  points  do  not  scatter  randomly  around  zero.  Rather, 
faint  sources  tend  to  be  brighter  in  the  ISO  7  pm  band,  whereas 
bright  sources  tend  to  be  brighter  in  the  IRAC4  band.  The  rela¬ 
tion  found  to  convert  between  the  two  bands  is 


mcigiscn  =  0.866  x  tnagjRA C4  +  1.024. 


(2) 


The  slopes  and  zero  points  of  this  linear  fit  are  similar  for  all 
fields,  with  a  cross-over  ( Spitzer  IRAC  4  equal  to  ISO  7  pm  mag¬ 
nitude)  between  670  and  7.n7.  The  trend  might  be  a  reflection  of 
a  dependence  of  the  strength  of  the  SiO  fundamental  band  at 
~  7.5  pm  and/or  a  water  band  at  ~  5.5  pm  on  the  brightness  of 
the  star.  At  least  the  SiO  first  overtone  band  at  ~  4.0  pm  has 
been  found  to  decrease  in  strength  from  semi-regular  variables 
to  Mira-like  variables  (Aringer  et  al.  1999).  On  top  of  that,  in¬ 
strumental  effects  in  one  or  both  missions  cannot  be  excluded. 


3.2.  Comparison  with  GLIMPSE-II 

The  Galactic  Legacy  Infrared  Mid-Plane  Survey  Extraordinaire 
(GLIMPSE-II;  Benjamin  et  al.  2003),  is  a  Spitzer  Legacy 
Science  Programme  of  most  of  the  inner  Galactic  disk  using  the 
IRAC  instrument.  Our  programme  has  overlaps  the  GLIMPSE- 
II  survey  in  the  fields  Bulge  1,  Bulge  2,  Bulge  4,  and  Bulge  N  1. 
A  comparison  with  this  large  data  set  is  especially  interesting  be¬ 
cause  the  same  instrument  and  the  same  filter  bands  have  been 
used,  though  the  observing  strategies  and  data  reduction  tech¬ 
niques  are  different. 

In  Fig.  [4]  we  show  the  magnitude  differences  between  our 
and  the  GLIMPSE-II  catalogue  for  all  four  IRAC  channels  for 
the  Bulge  N  1  field,  since  there  is  a  large  number  of  sources 


1.0 


is  _0-5  “ 

"o. 

1/1 

-i-or . ■ . ■ . ■ . ■ . ■  .  .  .  .1 

5  6  7  8  9  10 

Spitzer  IRAC4 

Fig.  3.  Magnitude  difference  between  the  Spitzer  IRAC  4  band 
and  the  ISO  7  pm  LW2  band  for  264  sources  in  common  in  the 
NGC  6522  field.  The  dotted  line  marks  the  identity  relation,  the 
dashed  line  is  a  linear  least  squares  fit  through  the  data  points, 
excluding  those  sources  fainter  than  973  in  either  of  the  data  sets 
to  avoid  Malmquist  bias. 


in  common  in  this  field  (31  191).  These  plots  look  very  similar 
for  the  other  fields  with  overlap.  In  general,  the  agreement  be¬ 
tween  our  and  the  GLIMPSE-II  catalogue  is  good  in  the  bright 
to  medium  brightness  range.  However,  at  the  faint  end,  start¬ 
ing  at  around  10m  -  1  lm,  we  notice  a  strong  trend  such  that  the 
magnitudes  become  increasingly  fainter  in  our  catalogue  than  in 
the  GLIMPSE-II  catalogue.  This  is  not  only  true  for  the  flag  B 
sources  (most  of  the  red  dots  at  faint  magnitudes  in  Fig.  [4|,  but 
also  for  flag  A  sources  fainter  than  this  magnitude  limit,  and  for 
all  filters.  The  same  trend  is  revealed  in  a  comparisorp] between 
the  GLIMPSE-II  catalogue  and  the  one  of|Ramfrez  et  al.  (2008) 
see  next  section).  According  to  B.  Babler  of  the  GLIMPSE-II 
team  (private  communication,  2008),  this  trend  at  the  faint  mag¬ 
nitude  limit  arises  because  of  GLIMPSE-II’s  limitation  to  “sin¬ 
gle  frame”  photometry.  GLIMPSE-II  does  photometry  only  on 
single  bed  frames  instead  of  mosaics,  because  each  patch  of  the 
sky  is  observed  only  twice.  This  strategy  may  cause  a  Malmquist 
bias,  such  that  GLIMPSE-II  magnitudes  will  be  increasingly  too 
bright  the  closer  a  source  is  to  the  faint  limit.  Thus,  for  sources 
at  the  faint  limit,  our  magnitudes  are  likely  more  accurate  than 
GLIMPSE-II  magnitudes.  Finally,  we  also  notice  some  satura¬ 
tion  effects  above  ~  675  in  IRAC  1,  and  a  few  sources  that  are 
significantly  fainter  in  the  GLIMPSE-II  catalogue  than  in  our 
catalogue  that  can  be  found  at  all  magnitudes. 

Due  to  the  large  number  of  sources  in  common  with  the 
GLIMPSE-II  survey,  we  are  able  to  check  whether  or  not  the 
combined  errors  are  of  realistic  magnitude.  To  do  so,  we  inspect 
the  sigma-factor,  which  is  also  discussed  in  the  aforementioned 
comparison  between  the  GLIMPSE-II  and|Ramfrez  et  al.  (2008) 
catalogues: 


our  mag  -  GLIMPSE  mag 

a  =  /  (3) 
y(our  error)2  +  (GLIMPSE  error)2 

If  the  combined  errors  are  close  to  the  magnitude  differences,  the 
sigma-factor  will  have  a  Gaussian  distribution  of  width  1.  If  the 
errors  are  underestimated,  the  distribution  will  be  broader  than 
that  and  narrower  if  the  error  estimates  are  too  large.  Figure  [5] 
displays  the  sigma  distribution  for  all  four  IRAC  channels  for  the 
field  Bulge  N  1,  along  with  a  Gaussian  fit  to  the  data.  Also  plot¬ 
ted  is  a  Gaussian  for  a  “perfect”  distribution  with  the  same  area 
under  the  graph,  i.e.  when  magnitude  differences  and  combined 
error  are  on  average  of  the  same  size  (cr  =  1).  The  source  selec¬ 
tion  was  restricted  to  the  magnitude  range  770-1070  to  avoid  the 
tail  of  faint  sources  with  a  Malmquist  bias  (see  above).  We  see 

1  http://www.astro.wisc.edu/sirtf/glm2_galcen_comparison.pdf 
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Fig.  4.  Magnitude  differences  between  our  and  the  GLIMPSE-II  catalogue  for  the  sources  in  common  in  the  N  1  field.  Sources  with 
quality  flag  A  in  our  catalogue  are  represented  by  black  dots,  all  other  quality  flags  by  red  dots. 


that  the  distributions  are  centred  on  zero,  but  there  is  a  surplus 
of  negative  sigma  values,  because  of  the  sources  with  negative 
magnitude  difference,  as  noted  above.  The  actual  distributions 
for  IRAC  1  and  2  are  only  slightly  broader  than  what  would 
be  expected  if  the  errors  were  correctly  estimated.  This  means 
that  the  combined  errors  for  these  channels  are  only  slightly  un¬ 
derestimated.  For  IRAC  3  and  4,  the  distributions  are  definitely 
broader  than  in  the  ideal  case.  For  these  channels  we  have  to 
assume  that  the  combined  errors  are  underestimated  by  factor  of 
1 .8  to  2.2.  Indeed,  our  error  estimates  are  on  average  smaller  than 
the  ones  estimated  by  GLIMPSE-II,  but  this  is  mostly  due  to  the 
longer  exposure  time  per  pixel  in  our  survey,  and  we  are  doing 
the  photometry  on  the  mosaicked  images.  If  we  used  two  times 
the  error  as  given  by  GLIMPSE-II  in  Eq.  [3]  the  sigma  distribu¬ 
tion  would  have  the  perfect  width  of  1.0  for  IRAC  1,  would  be 
slightly  too  narrow  for  IRAC  2  (too  large  errors),  and  would  still 
be  too  broad  for  IRAC  3  and  4  (too  small  errors).  We  therefore 
conclude  that  our  error  estimates  for  IRAC  1  and  2  are  proba¬ 
bly  only  slightly  too  small,  and  too  small  by  a  somewhat  larger 
amount  for  IRAC  3  and  4.  The  GLIMPSE-II  errors,  on  the  other 
hand,  are  probably  correctly  estimated  for  IRAC  1,  slightly  too 
large  for  IRAC  2,  and  slightly  underestimated  for  IRAC  3  and 
4.  Table  [6]  summarises  the  width,  centre,  and  peak  values  of  the 
Gaussian  fits  to  the  sigma  distributions  of  Fig.  [5] 


Table  6.  Width,  centre,  and  peak  values  of  the  sigma  distribu¬ 
tions. 


Band 

CT 

Centre 

Peak 

IRAC  1 

1.42 

-0.22 

241.5 

IRAC  2 

1.27 

0.17 

159.0 

IRAC  3 

1.80 

-0.24 

202.7 

IRAC  4 

2.15 

-0.17 

98.6 

3.3.  Comparison  with  Ramirez  et  al.  2008 ) 


The  overlap  with  the  GALCEN  survey  of  Ramirez  et  al.  (2008j> 
is  limited  to  the  Bulge  N  1  field.  The  reason  for  this  is  that 
the  Spitzer  time  allocation  committee  tried  to  avoid  redundant 
observations  of  the  same  fields  of  sky  within  GO  programmes 
(GLIMPSE-II,  on  the  other  hand,  is  a  legacy  survey).  A  compar¬ 
ison  between  our  own  and  the  GALCEN  observations  is  useful, 
because  for  our  N  2  field  we  will  have  to  use  the  IRAC  observa¬ 
tions  of  GALCEN.  Nevertheless,  because  IRAC  channels  1  and 
3  are  seeing  a  patch  of  the  sky  neighbouring  that  of  channels  2 
and  4,  we  got  a  small  overlap  with  Ramirez  et  al.|  ( 2008  i  in  our 
Bulge  N  1  field  in  order  to  fully  sample  that  field  in  all  bands. 
Adopting  a  stringent  170  search  radius,  there  are  5172  sources 
in  common  in  any  one  of  the  four  IRAC  bands.  In  the  Bulge  N  2 
field,  337  sources  detected  with  MIPS  24  have  a  counterpart  in 
Ramirez  et  al.  (2008),  adopting  a  less  stringent  370  search  radius 
because  of  the  less  precise  MIPS  coordinates. 
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RAC[3.6]  IRAC[4.5] 


Sigma  factor  (Prg  2345  -  GLIM PSE  —  II)  Sigma  factor  (Prg  2345  -  GLIM PSE  —  II) 


RAC[5.8] 


Sigma  factor  (Prg  2345  -  GLIM  PSE  —  II) 


RAC[8.0] 


Sigma  factor  (Prg  2345  -  GLIM  PSE  —  II) 


Fig.  5.  Distribution  of  sigma-factors  (Eq.  [3])  for  all  IRAC  channels  in  the  field  N  1  in  the  magnitude  range  7?’0  -  1  ()."().  The  solid  line 
is  a  Gaussian  fit  to  the  distribution.  The  dotted  line  indicates  a  Gaussian  of  width  1  and  the  same  area  under  the  graph,  the  ideal  case 
of  correctly  estimated  errors. 


Figure  [6]  shows  the  magnitude  differences  found  for  the 
sources  in  common  withjRarmrez  et  alT  (j2008[>  in  the  Bulge  N  1 
field.  The  agreement  in  the  IRAC  1  and  2  bands  is  very  good.  In 
the  IRAC  3  band,  the  GALCEN  magnitudes  (of  flag  A  sources  in 
the  medium  brightness  range  Tl’O  -  1  (CO)  are  on  average  ~  0'.n05 
fainter  than  our  magnitudes.  In  the  IRAC  4  band,  a  “knee”  ap¬ 
pears  for  sources  brighter  than  8™2.  Since  the  same  effect  is 
found  in  the  aforementioned  comparison  between  GLIMPSE- 
II  and  GALCEN  data,  the  origin  of  this  “knee”  most  probably 
comes  from  an  underestimation  of  fluxes  for  sources  brighter 
than  ~  8m2  in  IRAC  4  by  GALCEN.  We  will  have  to  con¬ 
sider  these  differences  between  our  IRAC  photometry  and  that 
of  Ramirez  et  al.  ( 2008|)  when  comparing  the  results  for  Bulge 
N  2  to  those  of  other  fields. 


The  differences  between  our  reduction  and  that  of  GALCEN, 
on  the  one  hand,  and  GLIMPSE-II,  on  the  other,  should  be  noted. 
While  in  our  reduction  the  background  is  determined  for  each 
pixel  in  a  45  x  45  pixel  box  centred  on  that  pixel  and  subtracted 
prior  to  the  photometric  measurements,  GALCEN  uses  the  local 
background  determined  by  PRF  fitting  and  subtracts  it  from  the 
corresponding  aperture  flux  of  each  detected  source.  For  the  fit¬ 
ting,  PRFs  pre-defined  by  the  SSC  are  used.  GALCEN  uses  the 
aperture  corrections  as  recommended  by  the  SSC,  but  adopts  the 
PRF  fluxes  for  their  final  catalogue  (except  for  sources  where 
the  ratio  between  PRF  flux  and  aperture  flux  exceeds  1.5;  for 
those  sources  the  aperture  flux  is  adopted).  We  adopted  the  PRF 


flux  in  our  catalogue,  derived  from  the  most  recent  mean  PRFs 
provided  by  the  SSC.  GLIMPSE-II,  on  the  other  hand,  performs 
photometry  on  individual  bed  frames,  not  on  mosaics.  A  com¬ 
bination  of  PRF  fitting  and  aperture  photometry  called  “tweak¬ 
ing^  is  used  to  measure  the  flux  of  a  source.  PRFs  dynam¬ 
ically  determined  from  each  individual  frame  are  fitted  to  the 
sources  and  subtracted  from  the  image.  On  the  residual  image, 
aperture  photometry  was  performed  around  the  positions  where 
PRFs  had  been  subtracted  in  the  previous  step.  If  the  resultant 
aperture  flux  is  substantially  positive  or  negative,  the  source  has 
been  under-  or  over-subtracted.  This  residual  aperture  flux  is 
then  subtracted  from  or  added  to  the  PRF  flux  to  compensate  for 
the  under-  or  over-subtraction.  However,  tweaking  was  applied 
rather  sparingly  and  only  in  the  IRAC  1  and  IRAC  2  bands  of  the 
GLIMPSE-II  observations.  These  differences  in  reduction  have 
to  be  kept  in  mind,  just  as  for  the  comparison  between  GALCEN 
and  GLIMPSE-II. 


3.4.  Comparison  of  MIPS  24  measurements  with  Hinz  et  ai. 


(2008) 

A  catalogue  of  Spitzer/MIPS  24pm  sources  towards  the  region 
of  ~  1:5  x  8°0  around  the  Galactic  centre  has  recently  been 
presented  by  Hinz  et  al.  (2008 1.  With  this  survey,  we  have  371 


2  http://www.astro.wisc.edu/sirtf/glimpse_photometry_vLO.pdf 
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Fig.  6.  Magnitude  differences  between  our  photometry  and  that  of 
as  in  Fig.  [4] 


Ramirez  et  al. 


2008  for  the  small  overlap  in  field  N  1 .  Symbols 


sources  in  common  (i.e.  unique  counterparts  within  a  one  phys¬ 
ical  pixel  =  2755  search  radius)  in  our  Bulge  4  field,  280  in  the 
Bulge  N  1  field,  and  129  in  the  Bulge  N2  field.  Thanks  to  the 
longer  exposure  time  of  our  observations  (33 1  s  per  pixel  com¬ 
pared  to  15  s  per  pixel  in  the  fast  scan  mode  used  by  Hinz  et  al.), 
our  catalogue  reaches  much  fainter  flux  levels. 

For  the  N  1  field,  we  only  have  a  small  overlap  with  the  sur¬ 
vey  of  Hinz  et  al.  (2008|>.  In  the  N  2  field  there  is  strong  diffuse 
galactic  emission,  and  a  small  fraction  of  the  sources  might  actu¬ 
ally  be  false  detections.  We  thus  concentrate  on  the  Bulge  4  field 
for  a  comparison  with  the  Hinz  et  al.  (2008|)  catalogue.  This  is 
also  the  field  with  the  most  sources  in  common.  In  Fig.  [7]  we 
show  the  magnitude  differences  between  our  measurements  and 
that  of  Hinz  et  ah](|2008]>  as  a  function  of  the  source  magnitude  in 
our  catalogue.  There  might  be  a  general  offset  of  ~  CP  15  in  the 
sense  that  our  magnitudes  are  on  average  fainter,  but  the  scatter 
is  considerable.  At  least  there  seems  to  be  no  trend  at  the  faint 
magnitude  end. 

As  for  the  comparison  with  GLIMPSE-II,  we  also  checked 
the  error  estimates  of  the  MIPS  24  data  by  investigating  the  dis¬ 
tribution  of  sigma-factors  (Eq.  [3}.  The  result  of  this  exercise  is 
that  the  combined  errors  are  too  small  by  a  factor  of  six  or  more. 


Even  adopting  twice  the  error  estimate  of  Hinz  et  al.  (2008), 


which  is  somewhat  larger  than  ours,  instead  of  a  combination 
of  the  errors,  gives  a  sigma  distribution  that  is  wider  by  a  fac¬ 
tor  of  four  than  what  is  expected  from  a  correct  estimate.  We 
thus  conclude  that  the  errors  are  still  largely  underestimated  in 


both  catalogues,  and  that  the  uncertainty  due  to  crowding  and 
strong  background  radiations  renders  a  precise  flux  determina¬ 
tion  in  the  vicinity  of  the  Galactic  centre  impossible.  Some  part 
of  the  found  differences,  however,  may  be  explained  by  a  real 
variability,  since  many  of  the  bright  sources  are  expected  to  be 
AGB  variables  in  the  bulge  (Glass  et  al.|2009]l.  At  a  wavelength 
of  20  fin l,  a  full  amplitude  of  up  to  077  is  found  for  C-rich  pul¬ 
sating  AGB  variables  in  the  solar  neighbourhood  (see  e.g.  Fig.  7 
of|Le  Bertre|1992]). 

4.  Colour-magnitude  diagrams  and  comparison 
with  theoretical  isochrones 

A  final  check  on  the  quality  of  our  data  reduction  are  colour- 
magnitude  diagrams  (CMDs)  in  combination  with  theoretical 
isochrones.  If  the  flux  measurement  has  no  systematic  error  (and 
also  the  isochrones  are  producing  realistic  colours),  the  observed 
giant  branch  will  be  described  well  by  the  isochrone.  This  check 
is  restricted  to  the  IRAC  bands,  since  for  many  sources  the 
MIPS  24  band  will  be  affected  by  dust  emission,  which  is  very 
difficult  to  model  accurately  in  isochrones.  We  chose  to  use  the 
most  recent  isochrones  from  Marigo  et  al.  (2008 1.  A  web-form 
can  be  used  for  computing  these  isochrones  with  different  pa¬ 
rameter^]  We  computed  an  isochrone  with  an  age  of  lOGyrs, 
metallicity  Z  =  0.019  (solar  metallicity),  and  no  dust  formation. 


http://stev.oapd.inaf.it/cgi-bin/cmd_2.! 
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4  6 

Our  magnitudes 


Fig.  7.  Magnitude  differences  between  our  MIPS  24  magnitudes 
and  Hinz  et  al.  ( 2008]),  for  our  Bulge  4  field. 


Table  7.  Coverage  of  our  fields  by  different  extinction  maps  of 
the  Galactic  centre  region. 


Field 

Ref. 

median  Av 

Bulge  1 

(1),  3 

4.38 

Bulge  2 

2,3 

2.45 

Bulge  3 

(1),  3 

2.00 

Bulge  4 

1,  3 

15.75 

Bulge  N  1 

1,  3 

7.25 

Bulge  N  2 

1,  3 

20.06 

NGC  6522 

1.  2,  3 

1.13 

References:  l:|Schultheis  et  al.|(|1999);  2:|Sumi|(|2004fr:  3:|Dutra  et  al,| 
(2003 1.  Numbers  in  brackets  mean  that  the  field  is  not  entirely  covered 
by  the  respective  map,  but  to  a  significant  fraction. 


In  the  mid-IR,  the  isochrones  are  not  very  sensitive  to  the  pre¬ 
cise  choice  of  the  parameters  such  as  age,  metallicity,  and  dust 
formation. 

To  disentangle  possible  systematic  errors  in  the  IRAC  chan¬ 
nels,  we  want  to  include  a  flux  measurement  that  is  independent 
of  our  reduction  method  in  the  comparison  with  the  isochrone. 
Since  all  our  fields  are  covered  by  2MASS,  which  can  be  re¬ 
garded  as  a  reliable  source  of  K$  magnitudes,  we  decided  to 
construct  dereddened  K$,o  vs.  (Ks  -  IRAC)o  CMDs. 

Most  of  the  fields  towards  the  Galactic  bulge  suffer  from 
strong  extinction.  Thus,  it  is  necessary  to  correct  the  photometric 
measurements  for  this  extinction  before  comparing  them  to  the 
theoretical  isochrone.  Several  extinction  maps  for  the  Galactic 
centre  region  have  been  published.  [Schultheis  et  al.  ( 1999)  de¬ 
rived  the  extinction  from  the  shift  of  the  giant  branch  in  the  Ks 
vs.  (J  -  K$)  CMD  based  on  DENIS  data.  The  extinction  values 
published  in  that  work  are  given  in  the  V  band.  Another  map 
of  Ay  values  towards  fields  in  the  Galactic  bulge  is  presented 
by|Sumi  (|2004|>.  These  are  based  on  the  extinction  measured  on 


red  clump  giants  in  fields  of  the  OGLE-II  survey  ( Udalski  et  al. 
2000  ).  Finally,  the  map  of  Dutra  et  al.  ( 2003 )  uses  the  same  tech¬ 


nique  as  Schultheis  et  al.  (|1999|  and  applies  it  to  2MASS  data. 
That  map  gives  the  reddening  value  in  the  Ks  band.  Its  spatial 
resolution  is  lower  than  that  of  the  other  maps,  but  it  is  the  only 
one  that  fully  covers  all  of  our  Spitzer  fields. 

Table[7]summarises  the  coverage  of  our  seven  fields  by  these 
extinction  maps,  and  the  median  extinction  in  the  V  band  is 
given.  The  NGC6522  field  is  the  only  one  that  is  covered  by 


all  three  extinction  maps.  Since  this  field  also  has  the  lowest  ex¬ 
tinction  of  all  our  Spitzer  fields  so  that  uncertainties  in  the  ex¬ 
tinction  determination  should  play  a  minor  role,  it  serves  here  as 
a  “fiducial”  field  for  the  comparison  with  theoretical  isochrones 
in  CMDs.  We  used  the  average  values  for  the  diffuse  inter-stellar 
medium  determined  by  Indebetouw  et  al.  ( 2005  i  for  the  extinc¬ 
tion  in  the  IRAC  bands:  A^s  :  A3 6  :  A45  :  A58  :  Ago  =  1  : 
0.56  :  0.43  :  0.43  :  0.43,  as  well  as  Av  :  A*, s  =  1  :  0.089. 
Using  the  somewhat  different  values  found  by  other  studies  (e.g. 
Flaherty  et  al.  2007 j  Nishiyama  et  al.|2009 1  has  only  a  negligible 


impact  on  the  results.  Figure  |8]  shows  CMDs  of  the  NGC6522 
field  involving  the  2MASS  Ks  magnitude  and  the  magnitudes  of 
the  four  IRAC  bands. 


The  following  observations  can  be  made  from  these  CMDs. 
At  the  faint  end,  the  distribution  in  colour  is  quite  broad  since  the 
errors  become  quite  large  for  these  faint  sources,  and  some  of  the 
objects  probably  are  background  galaxies.  There  are  also  a  small 
number  of  objects  with  extremely  red  or  extremely  blue  colours. 
Possibly,  false  identifications  are  involved  in  some  cases  (we  in¬ 
clude  only  sources  with  a  2MASS  counterpart  within  270  search 
radius  for  the  comparison  with  CMDs).  At  intermediate  bright¬ 
nesses,  the  distribution  in  colour  is  much  narrower,  and  only 
very  few  outliers  can  be  found.  At  the  bright  end,  the  number 
of  sources  far  from  the  isochrone  locus  increases  again.  This  is 
partly  due  to  saturated  sources,  e.g.  in  IRAC  1  (blue  “sequence” 
at  the  bright  end),  and  to  dust  emission  that  already  shows  up 
at  IRAC  4  wavelengths.  In  general  we  find  good  agreement  be¬ 
tween  the  observed  and  the  predicted  location  of  the  giant  branch 
in  the  CMDs.  In  the  (Ks  -  IRAC  l)o  CMD  (upper  left  panel  of 
Fig#  the  tip  of  the  giant  branch  in  the  isochrone  is  somewhat 
bluer  than  the  observed  tip.  In  the  (Ks  -  IRAC  3)o  CMD  (lower 
left  panel),  the  theoretical  isochrone  is  bluer  than  the  observed 
giant  branch  by  about  CPI  over  the  whole  brightness  range.  In 
this  filter,  the  agreement  with  the  GLIMPSE-II  survey  is  very 
good,  but  the  GALCEN  magnitudes  are  fainter  than  our  magni¬ 
tudes.  Thus,  by  adopting  the  GALCEN  magnitude  scale,  the  shift 
of  the  isochrone  would  be  somewhat  alleviated,  but  not  all  cata¬ 
logues  and  the  isochrone  can  be  made  agree.  We  do  not  find  any 
trends  in  magnitude  in  the  IRAC  1  filter  in  the  comparison  with 
any  of  the  other  surveys,  thus  we  assume  that  the  difference  is  re¬ 
lated  to  a  small  problem  with  the  isochrone  in  the  (DCs  -  IRAC  l)o 
colour.  The  same  might  be  suspected  for  the  (Ks  -  IRAC3)o 
colour. 


The  comparison  with  the  isochrones  also  allows  for  a  check 
on  how  deep  our  Spitzer  photometry  goes.  We  find  that  in  all 
IRAC  bands,  we  reach  sources  down  to  the  beginning  of  the 
He-core  burning  (horizontal  branch),  although  probably  with  a 
reduced  detection  probability  and  less  accurate  photometry.  The 
He-core  burning  phase  is  covered  until  its  end,  hence  also  the 
whole  AGB  evolution.  However,  the  sensitivity  of  the  present 
observations  is  still  far  too  low  to  reach  the  beginning  of  the 
RGB  phase,  and  our  investigations  will  be  limited  to  the  brighter 
half  of  the  RGB.  The  MIPS  photometry  does  not  reach  down  to 
the  horizontal  arm.  The  faintest  sources  that  are  detected  with 
MIPS  in  the  NGC6522  field  are  roughly  3m5  fainter  than  the 
RGB  tip  of  the  isochrones. 

As  for  the  colour  dependence  on  the  used  extinction  map, 
we  find  that  for  individual  sources  the  de-reddened  Ks  -  IRAC 
colour  may  vary  by  a  few  CP01 .  However,  no  general  trend  is  ob¬ 
served  when  different  extinction  maps  are  applied.  Thus,  though 
the  fine  structure  of  the  extinction  varies  somewhat  from  map 
to  map,  the  average  magnitude  of  the  extinction  is  very  similar 
among  the  maps. 
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Fig.  8.  Colour-magnitude  diagrams  of  the  field  NGC  6522  involving  the  2MASS  and  different  IRAC  magnitudes.  For  this  version, 
the  extinction  map  of  Schultheis  et  al.  ( 1999j)  has  been  used.  A  theoretical  isochrone  from  Marigo  et  al.  ( 2008)  is  included  as  a  blue 
line  (see  text  for  details).  Symbols  are  as  in  Fig.  [4] 


5.  Summary  and  conclusions 


We  present  a  catalogue  of  Spitzer  IRAC/MIPS  observations  of 
seven  fields  towards  the  Galactic  bulge,  sampling  a  range  of 
galactocentric  radii.  These  observations  allow  us,  amongst  other 
things,  to  study  the  mass  loss  of  a  large  and  homogeneous  sam¬ 
ple  of  RGB  and  AGB  stars  down  to  lower  luminosities  and  mass- 
loss  rates  than  previously  achieved. 


In  this  first  paper,  we  present  the  observations,  the  data  re¬ 
duction  procedure,  and  comparisons  to  other  mid-IR  surveys  of 
the  Galactic  bulge.  In  general,  we  find  good  agreement  with 
other  surveys.  The  comparison  between  the  Spitzer  IRAC  4  band 
and  the  ISOGAL  LW2  band  shows  good  agreement,  but  re¬ 
veals  a  slight  trend  with  magnitude.  GLIMPSE-II  magnitudes 
are  in  good  agreement  with  our  magnitudes  in  the  bright-to- 
medium  brightness  range,  but  a  strong  trend  is  present  at  the 
faint  end.  This  trend  is  probably  related  to  a  Malmquist  bias  in 
the  GLIMPSE-II  data  set.  The  error  estimates  of  GLIMPSE-II 
and  our  IRAC  photometry  are  reasonable  in  the  IRAC  1  and  2 
bands,  but  somewhat  too  small  in  the  IRAC  3  and  4  bands.  In  the 
comparison  with  GALCEN,  we  find  on  average  ~  O’." 05  brighter 
magnitudes  in  the  IRAC  3  band.  The  source  of  the  discrepancy  at 
the  bright  end  of  the  IRAC  4  band  is  probably  not  related  to  our 


catalogue.  A  comparison  with  the  MIPS  24  catalogue  of  Hinz 
et  al.  (|2Q08|>  reveals  that  our  magnitudes  are  probably  brighter 


at  the  level  of  ~  CP  15,  although  with  larger  scatter.  We  also 


find  good  agreement  between  our  data  and  recent  isochrones  in 
colour-magnitude  diagrams  for  at  least  three  of  the  four  IRAC 
bands.  The  ~  (P 1  offset  from  isochrones  involving  IRAC  3  de¬ 
serves  some  more  attention.  We  thus  may  assume  that  the  ob¬ 
servations  and  data  reduction  are  accurate  on  the  level  of  ~  (PI 
or  better,  as  well  as  precise  on  the  level  of  <  (PI,  except  for 
faint  sources  with  quality  grade  B.  The  science  exploitation  of 
the  data  will  follow  in  subsequent  papers. 
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Table  A.l.  Detected  and  extracted  sources. 


detected 

extracted 

SNR  <  5 

sources  outside 

both  criteria 

no  deblending 

37974 

37471 

380 

174 

51 

pass,  deblending 

37974 

37314 

483 

218 

41 

act.  deblending 

70474 

43655 

2145 

25885 

1211 

act.  &  pass,  deblend. 

49679 

39003 

1743 

9900 

967 

Appendix  A:  De-blending  or  not  de-blending? 


We  tested  the  comparability  of  the  fluxes  derived  by  PRF  fit¬ 
ting  and  those  resulting  from  aperture  photometry  for  the  IRAC 
3.6pm  channel  and  for  the  field  NGC6522  using  the  two  de¬ 
blending  capabilities  in  APEX.  The  first  method  is  passive  de¬ 
blending,  where  already  during  the  detection  process  multiple 
peaks  in  a  cluster  of  pixels  are  found  and  this  information  is  used 
in  the  PRF  fitting  module  to  simultaneously  fit  all  the  peaks.  The 
second  method  is  active  de -blending,  i.e.  the  number  of  assumed 
sources  is  iteratively  increased  during  PRF  fitting  to  decrease  the 
final  x2- 

We  use  the  point  source  probability  images  in  the  detection 
process  and  filter  the  raw  extraction  table  with  the  criterion  that 
SNR  is  greater  than  5  and  that  APEX  has  not  flagged  the  source 
as  a  point  source  outside  the  fitting  area.  This  excludes  very 
weak  sources,  sources  covered  by  only  one  image,  and  sources 
whose  position  is  not  correctly  identified.  The  numbers  of  de¬ 


tected  and  extracted  sources  for  the  runs  also  shown  in  Fig.  A.l 
are  given  in  Table  A.  1  along  with  the  number  of  sources  filtered 
by  the  two  criteria. 

The  active  de-blending  procedure  relies  on  the  x 2  values  to 
decide  whether  de-blending  is  necessary  or  not,  which  results 
in  heavy  de -blending  activity,  even  for  sources  whose  PRF  and 
aperture  fluxes  are  consistent.  We  therefore  disabled  active  de¬ 
blending  in  the  final  data  reduction. 

Passive  de-blending  uses  the  information  from  the  detection 
process  about  whether  different  entries  in  the  detection  table  be¬ 
long  to  the  same  cluster  of  pixels  above  the  detection  threshold. 
If  yes,  APEX  fits  the  cluster  simultaneously,  if  not,  the  single 
sources  one  after  the  other.  The  first  effect  of  switching  on  pas¬ 
sive  de-blending  is  that  the  number  of  faint  sources  that  lie  be¬ 
low  an  SNR  value  of  5,  hence  the  number  of  sources  rejected 
by  the  first  criterion,  increases  (see  Table  |AT|>.  Since  the  posi¬ 
tions  of  the  detected  sources  are  also  varied  to  fit  the  complete 
cluster,  the  number  of  sources  whose  pointing  is  outside  of  the 
tile  used  for  PRF  fitting  increases,  too  (Table  |A.1[).  The  num¬ 
ber  of  extracted  sources  is  thus  less  than  in  the  first  run  when 
passive  de-blending  is  not  activated.  Also  brighter  sources  are 
affected,  when  there  are  other  stars  within  5-10  pixels.  Then  the 
independent  stars  belong  to  the  same  cluster  and  are  fitted  simul¬ 
taneously,  although  the  PRF  of  one  star  does  not  contribute  at  all 
to  the  flux  in  the  small  tile  of  another  star  used  for  PRF  fitting. 
This  again  changes  the  position  of  the  centroid  of  the  central 
bright  star  in  the  cluster  and  also  the  flux  measured  by  PRF  fit¬ 
ting.  Since  the  displacements  are  smaller  than  the  pixel  size,  the 
aperture  flux  does  not  change  as  much  as  the  PRF  flux.  The  result 
is  that  the  bright  star  moves  away  from  the  line  of  good  agree¬ 


ment  in  the  fpRF-faperture-plane  (Fig.  A.l  i,  so  the  scatter  around 


this  line  is  increased.  The  same  effect  occurs  for  saturated  and 
super-saturated  sources  with  a  double  peak  leading  to  the  ob¬ 
served  spreading  of  sources  in  the  top-left  part  of  Fig.|AT|  Thus 
the  results  of  the  run  with  passive  de-blending  are  not  as  reliable 
as  the  results  of  the  run  without  any  de -blending. 
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Fig.  A.l.  Comparison  of  the  fluxes  resulting  from  PRF  fitting  and  aperture  photometry  without  de-blending  (first  row),  with  passive 
de-blending  (second  row),  with  active  de-blending  (third  row),  and  with  both  (last  row)  for  the  apertures  with  radii  of  2  pixels  (left) 
and  3  pixels  (right)  for  the  field  NGC6522  and  IRAC  1. 
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Appendix  B:  Magnitude  histograms  of  the  seven 
bulge  fields. 
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BLG1  field 


Fig.  B.l.  Histogram  of  magnitudes  for  the  Bulge  1  field.  The  bin  size  is  O'.'1 2. 


BLG2  field 


Fig.  B.2.  Histogram  of  magnitudes  for  the  Bulge  2  field.  The  bin  size  is  0."2. 


Fig.  B.3.  Histogram  of  magnitudes  for  the  Bulge  3  field.  The  bin  size  is  0."2. 
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BLG4  field 


Fig.  B.4.  Histogram  of  magnitudes  for  the  Bulge  4  field.  The  bin  size  is  O'.'1 2. 


Fig.  B.5.  Histogram  of  magnitudes  for  the  Bulge  N  1  field.  The  bin  size  is  022. 


NGC6522  field 


Fig.  B.6.  Histogram  of  magnitudes  for  the  NGC  6522  field.  The  bin  size  is  022. 
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Fig.  B.7.  Histogram  of  magnitudes  for  the  Bulge  N  2  field.  The  bin  size  is  CP2. 


